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In this article an efficient method to calculate the matrix elements of three-nucleon force is pre- 
sented. The new method is improved version of partial- wave decomposition of ref. [1], which sim- 
plifies expression to be evaluated and permits to reduce computational effort as much as six times. 
Proposed method naturally applies to Faddeev-type calculations but also can be used by any method 
relaying on partial-wave decomposition. 
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I. INTRODUCTION 

<N ; 

| Three nucleoli force (3NF) becomes an indispensable ingredient in realistic atomic nucleus calculations. In fact 
O , discussion is still vivid about the origin and structure of dominant three-nucleon force terms but even largest sceptics 
recognize its importance in nuclear physics [2]. Complexity of three-nucleon force turns to be a major obstacle towards 
its better understanding. Numerical implementation requires enormous analytical and programming effort. At the 
' same time even the best algorithms evaluating 3NF matrix elements become huge computer resources demanding tasks. 
In ref. [1] a nice Partial Wave Decomposition (PWD) technique of 3NF originating from two-pion exchange diagrams 
has been demonstrated. The present article improves this PWD by helping to simplify numerical implementation 
^3 \ procedure as well as reducing number of necessary floating point operations. 
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II. GENERAL REMARKS ABOUT PWD TECHNIQUE 

Three-nucleon force naturally decompose into three parts, which are identical up to a cyclic or anticyclic permutation 
' of the three interacting particles: 
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here in term particle k is considered as a formal spectator. While three-nucleon force Vijk is symmetric for 



exchange of any two nucleons, term is symmetric only for interchange of particles i and j. When working in 
isospin formalism,- since systems wave function is fully antisymmetric,- it is necessary to evaluate only one of three 
terms to get expectation value of the 3NF acting between three particles: 



<*|tS ifc |*)=3(* Vtf *). (2) 



Faddeev-type equations also require evaluation of single 3NF term, if they are written as follows: 

=G V l3 y + G V l f ) y, (3) 

(k) 

here is so called Faddeev component, Vij is two-nucleon interaction potential between nucleons i and j. Faddeev 
components are not fully antisymmetric, however they are antisymmetric for interchange of the nucleons in the active 
pair (i.e. component is antisymmetric for exchange of nucleons i and j). Fully antisymmetric systems wave 

function is a sum of three Faddeev components \& = + F^ + i*y . 

It is useful to work with Jacobi coordinates or momenta. One has three different Jacobi momenta sets, each selecting 
one of three particles as a spectator. I.e. if we select particle (k) as spectator, associated Jacobi momenta are defined 

by J^fc = \(k j — k i) and ~q \ = (| k k — 3 ), where fe,, kj and kk are individual momenta of the particles. This 
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momenta set is proper to Fadeev component . Following the standard procedure each Faddeev component 
is expanded in its proper bipolar harmonic basis \pqa) 
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here index a spans all the symmetry allowed combinations of the quantum numbers presented in the brackets: l p and 
l q are the partial angular momenta associated with respective Jacobi momenta; Sj and are the spins and isospins of 
the individual nucleons. Three nucleon system conserve total angular momentum J and its projection M. Individual 
isospins of the three nucleons couple to the total isospin T and its magnetic quantum number T z . 

Separation of the fully symmetric three-nucleon force into symmetric binary terms, as shown in eq.(l), can be done 

(k) 

in infinitely many different ways. In practice one associates the binary term with an expression obtained from 
the Feynman diagram where nucleon k is chosen as a formal spectator. Three-nucleon force related with two-pion 
exchange diagrams in its general form is usually written as [3]: 
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Cl ( ~T t . ~Tj) + C 3 {~T l . ~Tj) (~q l . Ifj) + C 4 (V fe . Q, X Qj^j (~T k . ~T t X ~Tj] 



where C\ , C3 and C4 are simply multiplicative constants and /j is scalar regularization function intended to cut out 
high- momenta contribution of the nucleon i; there is no strict regularization procedure to define cut-off function /, 
and several different regularization procedures have been used [3-5] . 
The last expression can be written in a compact form as: 



(6) 



with /( fe ) = (^.V,), I = (n.V.x^), = fifti.Qi), Wf - fifti.Q^Qu ^Wf = 

fj {j*o • Qj^j Qj x and fi = fi/ (^Qf + M^j. The momentum transfer Qi = fc ■ — k i is proper to Jacobi 

momenta set, where nucleon i is considered as spectator, giving: Qi = ~q' i — ~q i- 

In the following we will omit multiplicative constants C. Then all three terms of eq.(6) applied on the vector ^ 
regain the following form: 



q>' ~ W (i) IW (3)i S>. 

Huber et al. [1] proposed a nice method to evaluate these terms using the following expression: 
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here we have introduced convenient abbreviation J = 2 J + 1. In this expression J is the angular momentum of the 
three nucleon system and eventually of the bases a"", a'" and a. If potential term W is defined as vector (Wb and W 4 
cases) the angular momentum J' of the states a" and a' can be different from J. Last expression permits to evaluate 
matrix elements of the potentials W in their natural coordinate set, where Q is independent of one Jacobi momenta 
(p or q). This makes evaluation of matrix element i (p'q'a'\\ \\pqa) i quite rapid. Nevertheless one requires two 

/ \p'q'a') i , which consumes 
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and j {p"q"a"\ 



Jacobi basis transformation operations, namely k (p""q""a"" \p 
almost all the time required to estimate 3NF terms. Furthermore all the intermediate states (a' 



does not 



retain antisymmetry with respect to interchange of the nucleons in the active (coupled) pair. Therefore one is obliged 
to double the size of intermediate bases by including unphysical symmetric states, which are comfortably ignored 
when using Faddeev-type representation of systems wave-function. 
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III. IMPROVED PWD TECHNIQUE 

In this study we propose to redistribute 3NF terms in Faddeev equation, in the following manner: 



if ) = G ^* + -Go 



1 
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(9) 



Obviously these Faddeev equations satisfy the same three-particle Hamiltonian as eq.(3), thus they provide the same 
solution for the systems wave function '5. Nevertheless individual Faddeev components differ from components 
which satisfy eq.(3). 



Instead of using eq.(2), we calculate 3NF expectation value using identity: 



<*l^- fc |*) = ^(* 



Now eq.(7) is replaced with: 

2 L J 
Sum of first two terms (Ci and C3 ones) of 2-pion exchange force in eq.(6) we propose to evaluate as: 



(10) 
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Since systems wave function is completely antisymmetric, one has i {pqa — j {pqa . Furthermore W(Q j) is 
symmetric for exchange of the particles k and i, thus product j {p'q'a'\\ \\pqa)^ j {pqa \^) is a vector antisym- 
metric for exchange of the particles k and i. One can further demonstrate simple relation between matrix elements 
of k {p"q"a"\ JW) \p'q'a') i and k {p"q"a"\ J« \p'q'a') j : 



fc (p"g"a'V 0) IpV^ - R^'+Ch^^ k {p"q"a"\I^ \p'q'a') 1 , 



(13) 



where two phase factors reflect the symmetry of bra and ket states with respect to exchange of two active particles. Sec- 
ond phase factor is negative, since ket state is antisymmetric for (ki) pair. Furthermore matrices i {p'q'a' \ \pqa) { 
and j {p'q'a'\ \pqa) ■ are identical. Then it is easy to see that two terms in parentheses of cq.(12) cancel each 

other, unless the \p"q"a") k states are antisymmetric for exchange of the particles i and j (as should be, since vj k + V k f 

as well as arc symmetric in (ij) pair). When \p" q" a") k are physical, antisymmetric for exchange of the particles 
i and j, two terms in eq.(12) give identical contributions. Finally, one can write: 
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(14) 
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This expression contains only one particle-base transformation operation, which acts only between the physical states 
(ones being antisymmetric with respect to active particle pair). On the contrary old PWD eq.(8) required two particle- 
base transformation operations, which were coupling both unphysical (symmetric) and physical (antisymmetric) states. 
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Now we turn our attention to the C4 term of eq.(6), let evaluate it according to: 

/ III III /// |,T,'\ 
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x t (p'q'a'\\ U) Wf \\pqa), iipqa^) 

+ k ip"'q'"a'"\\ Wf ] \\p"q"a") k k ip"q"a"\ I \p'q'a') j 

X j ip'q'a'W^ Wf \\pqa) j j ipqa |*>) . (15) 

One should note the minus sign in front of the isospin term in the second line, which appears when splitting Vy fe 

into the product Iq(W^ .W W^). The matrix elements of the potentials (j ^W^ and are equal up to phase 

factor: 

i (p'?'a'|| (j) = {-)'*-" j ip'q'af* wf \\pqa) j (16) 

Term ( % )W^p is not anymore symmetric with respect to exchange of the particles k and i, therefore these matrices 
also couple physical (antisymmetric) states with the unphysical (symmetric) ones. Basically, phase factor in eq.(16) 
is positive for the transition between the physical states and is negative when one pass from physical state to an 
unphysical one. For the matrix elements of the isospin operators one has relation: 

k {p"q"a"\ Io \ P 'q'a') i = (-)*',+>',+*', k (p" q " a "\ I \p'q'a') j (17) 

Since sum of potential's terms vj k + V k p as well as term are symmetric with respect to exchange of the particles 
i and j the state \p" q" a") k present in the eq.(15) is antisymmetric for (ij) pair: thus first phase factor in eq.(17) is 
negative. Overall phase factor is negative if \p'q'a')- state is symmetric, whereas overall phase factor is positive if this 
state is antisymmetric with respect to pair (hi). Summarizing these remarks one can see that two additive terms in 
eq.(15) are identical, giving the matrix expression: 

k ip"'q"W" |*') = ^C- k ip"'q"W"\\ W[ k) \\p"q"a") k (18) 

x 
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k ip"q"a"\ I \p'q'a') j j ip'q'a'f ] wf \\pqa) . j ipqa |*> 

This time \p'q'a')- basis contain both physical and unphysical states, nevertheless other states {a, a" and a'") are 
physical. 

Nevertheless one can avoid arithmetics with unphysical states even evaluating term C4. In that aim we use identity 
~a k =25 — ~~o i — ~~o j with S representing total spin of three nucleon system to rewrite C4 term: 

C^wf . Wf y )l = C 4 (W^ . W~t ] )h + C 4 (W { b 3) . wf)I , (19) 

where wi^ = fj (lr j . Qj^j Qj x ^2 5 — 7?j^ and = fi (lr \ . Qij Qi x 7? j . The particular form of the last 

operator, as it is easily demonstrated in the Appendix, results in = 0. One can see that the term wi^ is 

symmetric for nucleon pair (hi). Therefore arguments used to derive eq.(14) can be followed to demonstrate: 

k ip"'q'"a"> |*') = - ]£' tlLl. ( fc (p"YW"\\ Wi k) \\ P "q"a") k (20) 

x k ip"q"a"\ Io \ P 'q'a') j j ip'q'a'\\ \\pqa). 3 ipqa |*>) 

This expression comprise only physical states and requires single basis transformation operation. Evaluation of 
3NF terms using old PWD of eq.(8) required two basis transformation operations: one j ip"q"a"\I |pV a ')i involving 
unphysical and physical states in both ket and bra parts (this transformation is numerically at least four times 
heavier than one coupling only physical states), whereas other k ip""q""a"" \p"'q"'a"') ■ transforming combination of 
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unphysical and physical ket states into the physical bra state (this transformation is numerically at least two times 
heavier than one coupling only physical states). As mentioned in [1] the basis transformation operations consumes 
essentially all the CPU time required to evaluate 3NF terms, while the time for the calculation of t (p'q'a'\\ \\pqa) i 
matrix elements is negligible. Therefore two-pion exchange 3NF terms can be evaluated by as much as six times faster 
using eq.(14) and eq.(20) than using old PWD. 

In the end I should remark that presented method can be applied to evaluate any momentum or coordinate space 
3NF, which has a form: Vijk = J2VijVjk- Application of similar PWD technique to evaluate UIX 3NF has been 

eye 

demonstrated in [6]. 

It is clear that described method is not directly applicable for more complex 3NF, which involves more than two 
meson exchanges, like Illinois force containing three-pion ring terms [8]. Nevertheless variation of this technique is 
possible to evaluate also the terms of type Vijk — VijVjkWki without introducing mixing with unphysical states. 
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APPENDIX A: 



Here we give explicit formulaes for the potential's terms W, required to evaluate two-pion exchange 3NF. When 
working with Faddeev equations one deals with binary particle interactions, which in proper Jacobi coordinate set 
are expressed as the functions of momenta p (and/or p') and which are independent of momenta q (q'). In contrary, 

our expressions for W^') -terms are the functions of Qj = kj — k j = ~~q'j — ~q j and they are diagonal in Jacobi 
momentum pj. In this sense evaluation of VF-terms are different from the evaluation of two-nucleon interaction terms, 
still they are similar in numerical complexity. 

The scalar potential term Wa^ = fj (jfj-Qj^ is the simplest case to evaluate: 



{p'q'a'\W^\pqa) j = ^^^,S jp ^ ~ 
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In the last expression Q = \/ q 2 + q' 2 — 2qq' u, P\(u) is Legendre polynomial. The f(Q) is model defined nucleon 
high- momentum cut-off function as introduced in eq.(5): in some models, as [4, 5], cut-off is set on momentum transfer 
Q = Q; in other models, as [3], cut-off can be set on single particle momenta Q = (q,q r ). One should note that the 
operator Wa^ conserves the three-nucleon angular momentum J, nevertheless it changes the parity of the state to its 
opposite. 

Next step is to evaluate reduced matrix elements of the vector potential = fj (~a j.~q j) ~q j: 
j(p'qW\\wl j) \\pq*)j = ^^Ki'^s^Kt'M' 2 ^ E ™ 

PP A k!+k 2 =K 




x(l + k 1 k 2 )Jp^J6k 1 k 2 T q V q jJ q JJ' { -f^ J '+^ 
k 2 A l q \ J j q j p J 1 J fc 2 lq A 

J \ J' 1 j' q J \ l'q h K 



q' kl q k2 9xK, (A3) 
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Finally, we present expressions for the reduced matrix elements of the vector potentials ^W^p = 
j (p'ct'a 



and W ( e 3) = fj (^j . Qj x (2S - ^)] : 
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and 
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The ( fc )Wjj" and ^potential terms conserve the parity of three-nucleon state, however it may change its 

total angular momentum. 

Using simple tensor algebra relations, one has: 

g h [ g h 1 J I- 9 > i 

In this equation the product {Tr j ® ^i}^ acts on spin one-half states and it is equal zero unless h = 0. However <? = 1, 
if ft = and |Qj (8) <5i| = as a vector product of two parallel vectors vanishes. Therefore the last equation is 

— > (i) ^ 

identically zero and Wy = 0. 

PWD basis transformation operations j (p"q"a"\ I \p'q'a') i have been discussed and necessary expressions have 
been given in rcf. [1, 7]. 
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